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The question of robusticity and the 
relationship between cranial size and 
shape in Homo sapiens 

Although cranial gracility is generally considered to be a characteristic feature 
of modern skulls, a quantification of the degree of development of cranial 
superstructures discloses varying levels of robusticity in certain recent and 
sub-recent populations. I n order to account for these structures, biomechanical 
interpretations in terms of masticatory stress and the effects of structural 
constraints have been put forward. This paper examines the multivariate 
relationships between robusticity and cranial dimensions in modern humans 
through canonical correlation and principal components analyses. These 
analyses reveal a very significant association between robusticity and cranial 
size, by which the larger the size of the skull, the greater the development of 
the cranial superstructures. In this context, size is biased towards maximum 
cranial length, various measures of upper facial breadth, basion-nasion length 
and palato-dental size. A statistically less important relationship between 
robusticity and cranial shape reveals that there are two main circumstances in 
which a modern skull will develop superstructures, in either very large crania, 
independent of the breadth of the vault, or in very narrow ones with large 
teeth. The cranial superstructures are not only strongly influenced by certain 
cranial dimensions, but are inter-correlated in their expression. These 
morphological features should not be treated as independent phylogenetic 
traits. 

The majority of recent crania are relatively small and gracile, clearly the 
result of a common temporal trend in the evolution of modern human cranial 
diversity. Although not unique in their level of robusticity among recent 
human groups, Australian aboriginals are the only population to show a 
relationship between robusticity and cranial breadth. It is suggested that the 
total morphological pattern observed in Australia does not reflect a plesio- 
morphic condition, but the specific history of occupation of that continent 
favouring size reduction with the maintenance of robusticity. 

Finally, it is shown that most early modern fossils are comparatively large 
and robust. It is suggested that the definition of a modern cranial form should 
be disassociated from the concept of gracility, only observed in some recent 
human populations. © 1996 Academic Press Limited 
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Introduction 

The term "robust" is used to describe something "strongly formed or constructed". This 
concept also applies to the human skeleton, which can be more or less robustly built. In the 
post-cranial skeleton, robusticity is normally identified in the overall size of the individual, the 
cortical thickness of long bones and the development of the areas of muscle attachment. 
Interpretations of the presence and development of post-cranial robusticity are generally 
associated with mechanical factors related to physical activity (Frost, 1967, 1972, 1982; Little, 
1973; Vaughn, 1981; Peck et al., 1984). In the cranium, robusticity is not only described by 
overall size and bone thickness, but it is also expressed by a number of distinct features, like 
ridges, tori and tubercles. These features have been called the "cranial superstructures" 
(Weidenreich, 1940), and their presence, absence and degree of expression in different 
populations of modern and fossil hominids has been the subject of many studies. Although 
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most of these studies offer only phylogenetic interpretations of these traits, there is a growing 
literature on studies of functional morphology, which includes works on the possible 
biomechanical origins of cranial superstructures. These works are always faced with the 
problems of identifying the ontogenetic and phylogenetic components within functional 
relationships, but they establish morphological units that co-vary in their growth parameters, 
and thus form units of phylogenetic change and inheritance. 

This paper intends to contribute to both interpretative aspects of cranial morphology- 
functional and phylogenetic— by further integrating cranial robusticity with metrical dimen- 
sions In order to do this, five aspects of the relationship between levels of robusticity and 
cranial size and shape will be examined— first, the character of such a relationship in terms of 
strength and constancy; second, the functional aspect of the variation observed; third, the 
regional distribution of cranial robusticity; fourth, the possible phylogenetic interpretations of 
both the relationship and the variation found; and last, implications for the current debate 
about modern human origins. 

Cranial robusticity in man 

M odern human crania vary markedly in size and shape. A large part of this variation can be 
attributed to populational differences, as groups differentiated in the course of their genetic 
response to environmental and demographic patterns through time. Accordingly, it is possible 
to use mathematical methods to identify crania from most modern populations on the basis of 
cranial measurements with a very high degree of certainty (H owells, 1973, 1989; Brace et al., 
1990; Pietrusewsky, 1990, 1992; Wright, 1992). H uman skulls vary also according to level of 
robusticity. As mentioned above, cranial robusticity is partly expressed in the presence and 
degree of development of a number of morphological features. Some of these features of 
robusticity are relatively common in modern crania, like well-formed supraorbital ridges, 
rounded zygomatic trigones, zygomaxillary tuberosities, occipital crests and occipital tori. 
Others are much more rarely observed, like rounding of the infero-lateral margin of the 
orbit, pronounced sagittal or frontal keeling, retromastoid processes, posterior supramastoid 
tubercles, mastoid tubercles and asterionic processes, and are sometimes found in specific 
regional groups. Still others, are typically observed in fossil hominid groupsand rarely or never 
in modern human skulls, like angular and supraorbital tori. 

Functional relationships 

As is the case in the post-cranial skeleton, the degree of robusticity of the human skull is also 
believed to be related to functional aspects of skeletal anatomy (H rdlicka, 1910; Washburn, 
1947; Moss & Young, 1960; Biegert, 1963; Aguirre & de la M acorra, 1964; Brace, 1964; 
Aguirre, 1968; Endo, 1970; Riesenfelt, 1974; O yen & Enlow, 1981; Russell, 1985; Shea, 1986; 
D ernes & Creel, 1988; H ylander et al., 1991). H owever, the precise nature of the functional 
complex that underlies cranial robusticity is still a matter of debate. 

M ost studies have concentrated on the mechanical role of the supraorbital ridges/ torus of 
hominids. Contrary to the occipital torus and occipital crest, the anatomical form of the 
supraorbital ridges/ torus is not directly defined by the insertion and/ or function of specific 
muscles or ligaments. T herefore, most functional interpretations consider this structure part of 
the stress resistance architecture of the skull. Endo (1965, 1966, 1970) demonstrated 
experimentally that the supraorbital area is mechanically deformed (bent) during biting. He 
postulated a model by which the facial skeleton acts as a rigid frame (beams) through which the 
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upward forces generated during biting are transferred to the supraorbital region, which also 
receives the downward bilateral pull of the temporal muscles. Endo further postulated that 
anterior dental loading was the main "stressing" force. Endo's studies were developed by 
Russell (1983, 1985), who found that it is the intermittent masticatory stress throughout an 
individual's life which stimulates cortical bone remodelling directed to acquiring a shape that 
counteracts the stress by minimizing the bending effect. Russell (1985) concluded that 
ontogenetic development of the supraorbital area will result in the bending stresses being 
neutralized by opposite structural curves. Russell did not include dental dimensions in her 
model, but D ernes & Creel (1988) showed that bite force depends on both the magnitude of 
the muscle force (which in turn depends on the specific dimensions that determine the load 
and lever arms of the contractile movement) and the occlusal pressure (dental area). T herefore, 
it is expected that dental dimensions will also be related to facial robusticity. Furthermore, 
Oyen et al. (1979a) found that rapid development of the supraorbital ridges coincides with 
dental eruption, and Oyen et al. (1979b) postulated that dental eruption changes the 
lever-arm :load-arm ratio, which changes the relation between masticatory muscles and 
the facial skeleton. H illoowala & Trent (1988a, b) developed thisfurther by arguing that it isthe 
force needed to elevate the mandible acting as a lever that determines the relationship between 
the supraorbital region of the skull and the masticatory apparatus through the pulling action 
of the anterior fibres of the temporalis. H owever, contrary to Endo and Russell, H illoowala & 
T rent (1988b) argue that it is not only anterior dental loading that creates the upward stressing 
forces, but that the posterior dentition is also involved. 

T here have been criticisms of E ndo's localized stress hypothesis (Picq & H ylander, 1989). A 
number of researchers have suggested that supraorbital tori develop not to counteract 
masticatory stress, but to fill the area between the brain and the face; in other words, 
supraorbital tori have a spatial role (M oss & Young, 1960; Biegert, 1963; Shea, 1986). This 
spatial model has received strong support from the morphometric studies of Ravosa (1988, 
1989, 1991) and the in vivo strain analysis in primates carried out by H ylander and co-workers 
(H ylander, 1984; H ylander et al., 1991). H ylander et al. (1991) find that the glabellar region of 
the frontal bone is indeed strained by mastication (equally by biting and chewing), but that this 
strain is very low, lower than in other portions of the facial skeleton. Although acknowledging 
that bone requires certain repeated levels of strain to complete its development, these authors 
do not accept that larger supraorbital tori are the result of structural adaptations to counteract 
greater masticatory stress. T hey argue that supraorbital ridges only perform a spatial function, 
a point that also explains the low levels of strain along the area. 

These different studies indicate that there is a functional relationship underlying supra- 
orbital development, whether directly associated to masticatory stress or to spatial arrange- 
ments, or both. However, many important points remain unresolved. The first is whether 
the features of robusticity like supraorbital and occipital tori, act as cranial superstructures 
with a common underlying relationship; second, whether the parameters underlying such 
a relationship are equally applicable in all modern populations; and third, whether it 
is ontogenetic or phylogenetic components of variation that determine most of the 
relationships observed. 

Temporal and spatial patterns 

Throughout the course of human evolution, features of robusticity like supraorbital and 
occipital tori, have been acquired, lost or changed in different groups, so that the recording of 
their presence, absence and character has been a useful tool in phylogenetic reconstructions. 
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As an example, it is possible to characterize populations of H omo erectus by the presence of a 
pronounced, horizontal supraorbital torus and supratoral sulcus on the frontal bone, while a 
morphological description of the Neanderthal population of the late Middle and Upper 
Pleistocene of E urope will certainly include the presence of a supraorbital torus that follows the 
contour of the orbits, arching over the middle portion of the supraorbital area. Modern 
humans, however, are not clearly characterized by their supraorbital morphology or by any 
other feature of robusticity as other hominid populations. Although most descriptions of 
modern H omo sapiens will include skeletal gracility as a characterizing trait, differing levels of 
robusticity are observed among recent and fossil modern skeletons. Early modern fossils have 
crania that are robust in relation to most recent populations— specimens like Omo I, Qafzeh 
9 and 6, Skhul 4 and 5, and the fragmentary remains from K lasies River M outh, show features 
like strong supraorbital ridges, which may form a real continuous torus, strong zygomaxillary 
tuberosities, developed zygomatic trigones and pronounced occipital tori. These traits have 
been mainly lost among recent populations, except in a few groups like Australian, Fueguian 
and Patagonian aborigines that continue to express high levels of cranial robusticity to the 
present. In these groups, cranial robusticity is interpreted as a plesiomorphic trait. Therefore, 
there has been a process of skeletal gracilization in modern humans that can be identified in 
most recent groups, although not a universal feature of the species, either temporally or 
geographically. These aspects of the evolution of a modern skull are intimately linked to the 
evolution of cranial diversity through the differential retention of robusticity in recent 
populations. 

Robusticity and cranial metrical dimensions 

As mentioned before, several studies have examined the relationship between supraorbital 
torus development and certain cranial dimensions (Russell, 1983, 1985; Demes& Creel, 1988; 
H illoowala & Trent, 1988a,b). In anatomical terms, the results of these investigations may be 
summarized as follows: there is an apparent relationship between dental size (anterior and 
posterior dentition), mandibular length, temporal fossae area and minimum frontal breadth 
with supraorbital torus/ ridge development. Furthermore, facial angle may also be correlated 
under conditions of great stress production (Russell, 1985). Lahr (1992, 1996) examined the 
relationship between the degree of expression of several features of robusticity, scored 
non-metrically, with metrical craniodental dimensions, and found that it is possible to predict 
the degree of expression of cranial robusticity in terms of cranial measurements with close 
to 100% accuracy. For many of the features under study, some of the cranial measure- 
ments responsible for most of the discrimination between the degrees of development of 
superstructures were the dimensions of the posterior dentition and of the cranial attachment 
of the temporalis, suggesting a close relationship between masticatory power and cranial 
robusticity. In spite of these results that would support the localized masticatory stress 
hypothesis, the presence of many other cranial variables in the discriminant functions 
highlights the role of cranial dimensions in the expression of cranial superstructures. 
Furthermore, like previous studies, Lahr (1992) analysed each of the features of robusticity 
independently, and although the metrical variables responsible for discriminating each trait 
were not exactly the same in each case, there were enough discriminating variables in common 
to suggest that at least the facial superstructures (supraorbital ridges, zygomatic trigone, depth 
of infraglabellar notch and zygomaxillary tuberosity) may respond to a single functional 
complex. 
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The following analyses attempt to address these issues by carrying out a more integrated 
study of the relationship between cranial measurements and robusticity. First, the existence 
of a close relationship between overall cranial robusticity (scored non-metrically for nine 
separate features) and craniodental metrical dimensions will be investigated, thus determin- 
ing the relative importance of the metrical and non-metrical traits studied in the morpho- 
logical variability observed. Second, the main underlying variables defining the relationship 
between robusticity and cranial dimensions will be examined for possible anatomical 
interpretations. Last, the geographical distribution of recent patterns of cranial robusticity 
will be described. 

M aterials and methods 

This study was carried out using a complete data set of cranial and dental measurements 
derived from several modern populations and a few individual fossils. M andibles were not used 
because they are too frequently missing, and the samples are unsexed. 0 nly for the discussion 
of the relationship between robusticity and sexual dimorphism, the levels of robusticity in 
males and females of two populations (Inuit and Fueguian) are compared. The number and 
composition of these two cranial series, representing restricted populations, enabled confident 
sexing of the skulls. T he exact number and composition of the samples used varied according 
to whether dental measurements were included or not, varying from 198 to 134 individuals in 
the first and second canonical correlation analyses, respectively. The principal component 
analysis (PC A) used the same sample as the second canonical correlation analysis. Variables 
with high frequencies of missing values were omitted, in order to maximize the number of 
crania in the samples. The geographical origin of these samples is described in Table 1. 

Two different sets of variables were obtained from each skull. The first one was composed 
of cranial measurements as described in Howells (1973), palatal length, a measure of the 
distance from the position of minimum frontal breadth to the face and the distance between 
left and right temporal muscles across the vault as described in Lahr (1992), and buccolingual 
dimensions of the first and second upper molars. I n order to maximize sample size, the dental 
dimensions were taken as the mean of left and right measurements, or either the left or right 
value when the tooth was present in one side only. This procedure is justified on the basis of 
the high correlation between left and right buccolingual values for the first and second molars 
(M 1: r = 09694, P <0-001, n = 182; M 2: r=09528, P <0-001) in the sample used (where teeth 
from both sides are present). T he second set of variables refers to the measurement of cranial 
robusticity, and it used ordinal grades to score degree of development. The features of 
robusticity examined were the development of the supraorbital ridges, zygomatic trigone, 
zygomaxillary tuberosity and sagittal keeling, the depth of the infraglabellar notch, degree of 
rounding of the infero- lateral margin of the orbit, and categories of nasal saddle, occipital 
torus development and occipital crest (Figure 1). T he grades and scoring procedure of each of 
these traits may be found in Lahr (1992, 1994, 1996). The complete set of metrical and 
non-metrical variables used is listed in Table 2. 

Statistical analyses 

The study of the relationship between a set of metric variables, representing cranial 
dimensions, and a set of non-metric variables, representing robusticity, was carried out 
through canonical correlation analyses. Canonical correlation methods will reveal the 
correlation underlying metric and non-metric variables as two sets, and answer the question of 
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Table 1 Composition and size of the samples used in the canonical 

correlation (CC) and principal components (PC) analyses 



R egion 


1st CCA 


2nd CCA/ PCA 


Europe 


19 


13 


North Africa 


7 


3 


Africa south of the Sahara: 






East Africa 


6 


7 


West Africa 


9 


7 


South Africa (bushman) 


6 


5 


India 


2 


1 


Southeast Asia 






Southeast Asia 


20 


14 


Andaman Islands 


4 


3 


East and Northeast Asia 


25 


12 


Inuit 


5 


4 


American Indian: 






North American Indian 


1 


1 


South American Indian 


3 


— 


Fueguian/ Patagonian 


25 


12 


Polynesian 


6 


3 


Island M elanesia 


5 


2 


Australia 


32 


27 


Mediterranean Epi-Palaeolithic: 






Afalou 


13 


9 


Taforalt 


3 


4 


N atufian 


3 


3 


Pleistocene fossils: 






Zhoukoudian UC (China)* 


2 


1 (U C 103) 


Liujiang*t 


1 


1 


Ohalo (Israel) 


1 


1 


Nahal Ein Gev (Israel) 


1 


1 



*M easured from casts at The Natural H istory M useum, London. 
tThe Liujiang fossil has no precise date (Brown, 1992). 



how the sets can be combined to make the correlation between components of the two sets a 
maximum. If the first canonical correlation is statistically significant, there is at least one way 
in which the domains of metric and non-metric variables are related (Cooley & Lohnes, 1962). 
Most users will want to examine the relationship between the original variables and the 
canonical vectors, but as Wherry (1984: pp. 300-301) points out, some workers have caused 
confusion by trying to interpret the sets of variable weights to interpret the results. The correct 
procedure is the use of the structure coefficients, which are the loadings of the variables on the 
canonical vectors of the sets. 

Redundancy coefficients of Stewart and Love are an alternative measure of association 
between two sets of variables. U nlikethe canonical correlation coefficients, these measures are 
asymmetrical with respect to the two sets. Redundancy coefficients measure the proportion of 
variance of the variables in the relevant set that is accounted for by the corresponding 
canonical variate. As Tatsuoka (1988) points out, redundancy coefficients can play an 
important role in modifying enthusiasm about the canonical correlation coefficient. We can 
think of redundancy coefficients as measuring the overlap between the two sets of variables 
(Timm, 1975). The redundancy measure can also be interpreted as a measure of the extent to 
which one set of variables (either metric or non-metric) can be reconstructed from knowledge 
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ZT-Zygomaxillary 
tuberosity 




(in c is u r a 
m alaris) 



(Malar tubercle) 



RO-Rounding of 
the infero-lateral 
margin of the orbit 



Figure 1. Cranial superstructures. In bold, the features of robusticity scored as ordinal grades and used in the 
present study. In parentheses, other tubercles and tori of the human skull not included in the present analysis. 



about the other set (G reen, 1978). A textbook by G reen (1978) provides a clear and integrated 
description of canonical correlation analysis, with an example worked from the raw data (not 
just from the correlation matrix). 

T he much better known technique of PC A (G reen, 1978) was performed on the same data 
set and using the same variables as the second canonical correlation analysis. PC A ignores the 
differences in variable types (whether metric or nonmetric) and seeks only to maximize 
explained variance, regardless of the a priori sets of variables. Given the different aims of 
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Table 2 Metrical and non-metrical variables used in the study 



M etrical 

(a) Description in H owells (1973) 



GO L— maximum cranial length 
BNL— basion-nasion length 
XFB— max. frontal breadth 
ASB— asterionic breadth 
NPH — nasion-porion height 
N LB- nasal breadth 
OBB- orbital breadth 
M A B— maxillary breadth 
BPL— basion-porion length 
SSS— subspinale subtense 
NAS— nasion subtense 
FRC — frontal chord 
PAC — parietal chord 
OCC — occipital chord 
GLS— glabellar subtense 
FRF— frontal subtense-fraction 



NOL— naso-occipital length 
BBH — basi-bregmatic height 
XPB— max, parietal breadth 
AUB— biauricular breadth 
NLH- nasal height 
OBH- orbital height 
JUB-bijugal breadth 
WM H — min. malar height 
ZM B— bimaxi Mary breadth 
FM B-bifrontal breadth 
EKB-biorbital breadth 
FR S — frontal subtense 
PAS— parietal subtense 
OCS— occipital subtense 
SOS— supraorbital subtense 
OCF— occipital subtense-fraction 



PA F— parietal subtense-fraction 
(b) Description in Lahr (1992) 

WIOB: minimum interorbital breadth 
PALL: palatal length 

WFBZFS: distance position of minimum frontal breadth and the zygomatico- frontal suture 
SI L L PA : distance between left and right temporal muscles across the vault 
BLM 1: buccolingual upper first molar length 
BLM 2: buccolingual upper second molar length 
*Non-metrical (features of robusticity) 

ST: development of the supraorbital ridges and formation of a supraorbital torus 

T R : development of the zygomatic trigone 

IN: depth of the infraglabellar notch 

N S: categories of nasal saddle 

ZT: development of a zygomaxillary tuberosity 

RO: rounding of the infero- lateral margin of the orbit 

SK : degree of sagittal keeling 

OT: configuration of the superior and supreme nuchal lines and formation of an occipital torus 
OCR: development of an occipital crest 



*D escription of grades and scoring techniques in Lahr (1992, 1994, 1996). 



canonical correlation and PCA, a comparison of the results obtained through these two 
techniques provides a good test of the findings. 

The program for canonical correlation analysis is written by one of us (R.S.V.W.), as an 
auxiliary program for the M V-N UTSH ELL package for archaeological analysis. The program 
for PCA comes from the distributed version of M V-N UTSH ELL (Baxter, 1994; Wright, 
1994). 



T he relationship between cranial robusticity and cranial dimensions 

T he results of the first canonical correlation analysis on 198 crania using all the non-metrical 
features of robusticity and a set of cranial metrics show that the non-metrical and metrical 
variable sets are correlated: canonical correlation of factor 1:0-85. This correlation is highly 
significant (% 2 = 525'74; P <0 0001). The variable loadings show that the first canonical vector 



Results 
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Figure 2. Plot of variables along the metric and non-metric vectors of the first canonical variate correlating 
cranial measurements and overall robusticity (measured as grades of expression of nine morphological 
features). Area in dotted line, metrical and non-metrical variables moderately to highly correlated with the 
first canonical variate. 



relates overall cranial size with intensity of non-metrical traits (robusticity), i.e. the bigger the 
crania the more robust and vice-versa (Figure 2). The metrical traits that are moderately to 
strongly correlated with the first canonical vector (correlations larger than 0-4) are maximum 
length (GOL, NO L), upper facial breadth (FMBJUB, EK B, 0 BB), facial-basal length (BNL), 
palatal size and protrusion (M AB, BPL), and occipital flexion (OCS). The variable loadings for 
the features of robusticity vary from 0-203 in the case of categories of nasal saddle to 0732 in 
the case of development of the supraorbital ridges, i.e. there is a lot of variation in the level of 
correlation of the various non-metrical features used. The loadings on the first canonical 
variate are monopolar, i.e. they are either positive or close to zero (those close to zero are 
not correlated with the first vector). These results show that the main correlation of cranial 
robusticity is not with a metrical shape factor, but with cranial size (chiefly represented 
by length and upper facial breadth). The redundancy values show that 24-9% of the 
non-metrical variance is explained by the metrical vector, whereas only 9-8% of the variance 
in metrical dimensions is accounted for by the non-metric vector. These values indicate that 
the direction of influence is from the metrics to the non-metrics, i.e. robusticity increases 
as cranial size increases (length and upper facial breadth), but the reverse effect is much 
weaker. 

The discussion now turns to interpreting the crania plotted on the scattergram. For the 
purpose of exposition, somewhat arbitrary lines are drawn to define "groups". These groups 
are examplars, not groups derived by any formal technique such as cluster analysis. 

Two anatomical groups defined by the relationship between cranial size and robusticity can 
be identified, one group composed of large and robust crania, and another one of small and 
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gracile skulls, circled in the upper right and lower left quadrants of Figure 3 respectively. It is 
also possible to identify two groups of crania that depart in opposite ways from the 
relationship, i.e. crania that are either too robust or too gracile for their size. As shown by the 
mean value of cranial and frontal breadth in both the robust and gracile groups (Graciles: 
XPB: x = 137 mm, XFB: x = 115 mm; Robusts: XPB: x=136 mm, XFB: x=114 mm), cranial 
breadth does not play any role in the definition of the observed relationship between 
robusticity and size. 

The relationship between facial robusticity and craniodental dimensions 

In order to examine whether there is a closer functional relationship underlying the 
development of only the facial superstructures (excluding features like occipital torus and 
occipital crest from the non-metric set of data), a second canonical correlation was performed. 
This included only features of facial (development of supraorbital ridges, zygomatic trigone, 
zygomaxillary tuberosity, rounding of the infero-lateral margin of the orbit and depth of 
infraglabellar notch) and vault (sagittal keeling) robusticity. Because of the potentially 
important role of masticatory stress in the degree of development of these features, two 
measures of dental size (buccolingual M 1 and M 2), palatal length, a measure of proximity of 
left and right temporal muscles across the vault and a measure of the position of minimum 
frontal breadth relative to the face were included in the metrical set of variables. The inclusion 
of dental measurements in the analysis unavoidably reduced the sample size. T his analysis was 
performed on a data set of 134 skulls. 

The canonical correlation between cranial metrical dimensions and facial robusticity was 
very strong (canonical correlation of factor 1:0-90). All the features of robusticity examined 
(ST, TR, IN, ZT and SK) show moderate to strong positive loadings on the first canonical 
variate, suggesting that they follow a common underlying relationship with the metrical 
variables. T he metrical variables with moderate to high loadings on the first vector are similar 
to the ones found in the first analysis [Figure 4(a)]— maximum length (GOL, NOL), maxillary 
breadth (MAB), upper facial breadth (OBB, FM B, EKBJUB), facial-basal length (BNL), but 
also include specific features of the supraorbital region, like glabellar protrusion (GLS)and the 
distance between position of minimum frontal breadth and the zygomatico- frontal suture 
(WFBZFS), and palatodental dimensions (PALL, BLM1, BLM2). The redundancy results 
indicate that the relationship obtained between the metrical dimensions and the facial features 
of robusticity is stronger, 37-2% of the variance in robusticity is explained by the metric set of 
variables, and 11% of metrical variation is explained by variation in facial robusticity. 

The second canonical correlation has a value of 077, and it is also significant (% 2 = 34072; 
P <0-0001). This second canonical variate is orthogonal to the first, and naturally describes 
crania that have a very different relationship between facial robusticity and metrical 
dimensions. The variables with moderate positive loadings (0'3-0-6) on the first factor of the 
second canonical variate are only two: sagittal keeling (SK ) and orbital height (0 BH ), and with 
negative loadings ( — 0-3- -0'6) are minimum interorbital breadth (W 1 0 B), distance between 
left and right temporal muscles across the vault (SILLPA)and rounding of the infero-lateral 
margin of the orbit (RO). The variables correlated significantly on the second factor of the 
second canonical variate are the same as on the first, with the difference that the loading of 
sagittal keeling is higher. The correlation of variables along the second canonical variate is 
bipolar, i.e. there are both positive and negative high loadings, and indicates that this variate 
identified a "shape" factor associating keeling of the vault and certain dimensions [Figure 
4(b)]. The metric factor accounts for 7'3% of the variance in the non-metric variables, while 
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Figure 3. (a) Plot of crania along the first canonical variate, separating robust crania on the upper right-hand 
quadrant from gracile ones on the lower-left quadrant. H istograms and box metrical values refer to the 
crania circled on the two extremes of the plot (the robust and gracile groups). The histograms [(b) gracile 
crania; (c) robust crania] show the expression of robusticity of the six non-metric grade features in these two 
groups of crania, ranging from very pronounced (black) to absent (white). The boxes show the mean values 
of the traits moderately to strongly correlated with the metric vector, whereas the mean values of X PB and 
XFB illustrate the absence of correlation of these variables with the first canonical vector obtained. 
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Figure 4. (a) Plot of variables along the metric and non-metric vectors of the first canonical variate correlating 
cranial-dental measurements and craniofacial robusticity (measured as grades of expression of six 
morphological features). Area in dotted line, metrical and non-metrical variables moderately to highly 
correlated with the first canonical variate; (b) plot of variables along the metric and non-metric vectors of the 
second canonical variate correlating cranial-dental measurements and craniofacial robusticity (measured as 
grades of expression of six morphological features). 
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the non-metric factor accounts for 1-8% of the variance in the metric variables. Therefore, this 
second canonical variate identified a very specific relationship that discloses a cranial shape 
associated with the development of keeling of the vault: pronounced sagittal keeling is 
correlated significantly with high orbits, narrow inter-orbital widths, spatially close temporal 
muscles across the vault and rounding of the infero- lateral margin of the orbit. Taken 
together, the relationships identified in the first and second canonical variates account for 
44'5% of the variation in craniofacial robusticity. 

As with the first canonical correlation analysis, the results of this second analysis using 
only facial superstructures in the non-metric vector allow the identification of facially robust 
and gracile skulls (Figure 5). The individual composition of these two groups is very similar 
to the robust and gracile groups of Figure 3, although not identical. Both groups (robust 
and gracile) of Figure 5 contain individuals that were outside the circled areas in Figure 3, 
but on the other hand, except for a single skull, all the crania that were classified as either 
robust or gracile in Figure 3 and used in the second canonical correlation analysis, are 
within the facially robust or gracile groups of Figure 5. This would suggest that most skulls 
show a similarly distinct pattern of cranial and facial robusticity (either gracile or robust), 
but that in a few skulls facial robusticity or the absence of, is more pronounced than the 
overall cranial pattern. This further reinforces the view that the degree of expression of 
each feature is correlated with the individual's level of robusticity, although certain groups 
of features like the facial traits analysed in the second canonical correlation, form a 
recognizable morphological complex. 

PCA analysis 

A PCA was performed using the same data and variable set as the second canonical 
correlation. The results of this analysis are very similar to those obtained in the canonical 
correlation analyses. The metrical variables that have a strong positive loading on the first 
principal component (GOL.JUB, FMB, EKB, NOL, MAB, BNL, OBB, ZMB, AUB, BBH, 
FRC) are mostly the same that have strong or moderate loadings on the first canonical variates 
of the canonical correlations, although the features of facial robusticity and sagittal keeling are 
only moderately correlated with the first principal component. This component accounts for 
27-3% of the total variance, and may be said to represent cranial size as defined before in 
terms of maximum cranial length, upper facial breadth, palatal size and facial-basal length, 
and a few other main cranial dimensions, like basal breadth, cranial height and frontal length. 
As was the case for the first canonical variate, the loadings on the first principal component are 
mainly positive or close to zero, confirming that the main structure in the data is one of metric 
size correlating with robusticity and the reverse. 

The second principal component, which accounts for 11-4% of the variance, is different in 
character from the first in the sense that the variables have both positive and negative loadings 
on this factor (like the second canonical variate). I n other words, there is a shape component 
involved. The variables with moderate to strong positive loadings on PC 2 are all related to 
width dimensions (X PB,XFB,SILLPA,AUB), whereas those with moderate negative loadings 
are related to robusticity, palatodental size and protrusion, and position of minimum frontal 
breadth (TR, IN, RO, ST, BLM 1, BLM2, SSS, PALL, WFBZFS). Therefore, the "shape" 
relationship disclosed by the second PC is different from that identified in the second canonical 
variate. It contrasts broad and gracile crania with robust and narrow ones, and it shows the 
close association between degree of facial robusticity and dental size. Orthogonal plots of the 
variables along principal components 1 and 2 show the features described above (Figure 6). 
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Figure 5. (a) Plot of crania along the first canonical variate, separating facially robust crania on the upper 
right-hand quadrant from facially gracile ones on the lower-left quadrant. H istograms and box metrical 
values refer to the crania circled on the two extremes of the plot [(b) gracile and (c) robust groups]. The 
histograms show the expression of robusticity of the six non-metric grade features in these two groups of 
crania, ranging from very pronounced (black) to absent (white). The boxes show the mean values of the traits 
moderately to strongly correlated with the metric vector, whereas the mean values of X PB and X FB illustrate 
the absence of correlation of these variables with the first canonical vector obtained. 
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The anatomical interpretation of this plot is complicated by the correlation of robusticity on 
both the first and second principal components, so that it is possible to find both broad and 
narrow skulls that are relatively robust (Figure 7). 

In order to interpret the spatial distribution of skulls along the first and second principal 
components, a combination of the size (represented by cranial length, upper facial breadth and 
facial-basal length) and shape variables was used to delineate arbitrary groups of crania with 
corresponding metrical dimensions (Figure 7). Table 3 shows the mean values of the metrical 
variables used to define these groups (those strongly correlated with PC 1 and PC 2) for the 12 
groups obtained. The most peripheral groups present comparatively extreme morphologies, 
whereas the four central groups (C , E , H , K ) are less differentiated. T he four internal groups 
comprise 64'9% of the crania studied. Figure 8 shows the gradient of robusticity visible from 
the upper left quadrant of Figure 7 (groups A and D ) towards the lower right one (groups G 
and I ). 1 1 is clear that the broad short skulls show the smallest development of facial robusticity, 
whereas both the very narrow, but not particularly large, skulls of group G and the long and 
narrow skulls of group I, show a very pronounced expression of facial robusticity. 

R egionality in recent cranial robusticity. The regional distribution of overall cranial robusticity can 
be examined if the crania in Figures 3 and 5 are identified by population of origin. T he group 
of skulls circled in the upper right quadrant of Figure 3, i.e. robust crania, are 23 in all: 15 
Australians (65-2%), three Fueguians/ Patagonians (13-0%), two Afalou/ T aforalt (8-7%), one 
Ainu (4'3%), one Maori (43%) and Zkn UC 101 (43%). The regional composition of the 
group of skulls with pronounced facial robusticity (circled in the upper right quadrant of Figure 
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5) isvery similar, 65-2% Australians, 21-7% Afalou/ T aforalt, 8-7% Fueguian/ Patagonian and 
4'3% Ainu. Among the crania along the axis of the relationship, which are neither particularly 
robust or gracile, those towards the upper right quadrant are mainly Australian, Fueguian/ 
Patagonian or Afalou/ T aforalt. These populations can be said to be comparatively robust 
among recent modern humans (Figure 9). The group of skulls circled in the lower left quadrant 
of Figure 3, i.e. small and gracile crania, are 44 in all: 13 Africans (29-5%), 13 Eastern Asians 
(29-5%), 11 Southeast Asians (25-0%) and seven Europeans (15-9%). The crania identified as 
gracile in terms of facial features (circled in the lower left quadrant of Figure 5) are 34 in total: 
11 African (32-3%), nine Southeast Asian (26-5%), six East Asian (17-6%), three Europeans 
(8'8%), two Inuit (5'9%), one Indian (2'9%), one Australian (from the northwest- 2'9%) and 
the fossil of Nahal Ein Gev,from Israel. These populations, the majority of recent populations, 
may be considered relatively gracile within present day cranial variation. 

The regional origin of the crania identified as either too robust or too gracile for their 
size in Figure 3, i.e. crania that do not follow closely the relationships size/ robusticity 
identified in the canonical correlation, does not help to understand why they depart from 
the main pattern. Among those skulls that show greater robusticity than expected for their 
size, we find two Southeast Asians, one African, one Amerindian, four Europeans, one 
Natufian, one Fueguian/ Patagonian, and three Australian. Among those that may be said 
to be relatively large for their gracility we find three Afalou/ T aforalt, three Fueguian/ 
Patagonian, Zkn UC 103, four Africans, and one Polynesian skull. The fact that some of 
the individuals that are too robust or too gracile for their size in Figure 3 fall within the 
robust or gracile groups of facial features in Figure 5, would suggest that certain individuals 
show different patterns of gracility/ robusticity for different functional complexes, and thus 
depart from an overall cranial pattern. 

Sexual dimorphism 

Sexual dimorphism is an important factor in the expression of robusticity, and it may affect the 
interpretations presented here in two ways. First, it is possible that the relationship between 
cranial dimensions and robusticity varies between males and females. T he crania used in this 
study are unsexed, and therefore, this possibility cannot be assessed. Nevertheless, the 
relatively small dispersion of data points, for each regional group, along the canonical 
correlations would suggest that the male and female crania contained in the sample behave in 
a similar way. 0 n the other hand, it is possible that biased regional samples towards males or 
females are affecting the identification of populations as robust or gracile, while if the regional 
samples were sexually balanced maybe each population would cover the range of variation in 
robusticity. In order to investigate this possibility, two sexed series, one a relatively gracile 
group, I nuits, and the other a robust population, Fueguians, were examined. T he skulls were 
sexed using features of the mandible, orbital sharpness and whenever possible, post-cranial 
information. Each skull was given a "robusticity score" representing the sum of the grade 
scores of development of the supraorbital ridges, zygomatic trigone, infraglabellar notch, 
zygomaxillary tuberosity, sagittal keeling, frontal keeling, occipital torus, occipital crest and 
external occipital protuberance, and the variation in the two populations compared (Figure 
10). From these data it is clear that females are more gracile than the males of their own 
population, although there is overlap between the sexes (more in I nuits). H owever, the gracile 
and robust character of each group is maintained. M ale I nuit skulls show "robusticity scores" 
similar to female Fueguian crania and no Inuit skull is as robust as Fueguian males. Plotting 
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Figure 10. Differences in robusticity between male and female Eskimo and Fueguian skulls, and comparison 
to the distribution of levels of robusticity within modern populations (Europe, sub-Saharan Africa, Southeast 
Asia, East Asia, Australia, Afalou-Taforalt, Natufian) and certain Upper Pleistocene modern fossil specimens. 
Robusticity expressed as an individual score representing the sum of the ordinal grades of development of the 
supraorbital ridges, zygomatic trigone, infraglabellar notch, zygomaxillary tuberosity, sagittal keeling, frontal 
keeling, occipital torus, occipital crest and external occipital protuberance. 
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the "robusticity scores" of the unsexed populations used in this study (the complete 
series) further confirms the identification of gracile (sub-Saharan African, Southeast 
Asian, East Asian, European, Natufian) and robust (Australian, and to a less extent 
Afalou/Taforalt) groups, with fossil European crania falling in the robust end of the scale 
(Figure 10). 

Figure 11 shows the identification by population of the crania along the orthogonal plot 
of the first and second principal components (Figure 7). H orizontally, along the axis of PC 1 
(correlated with size— maximum cranial length: GOL; bijugal, bifrontal, biorbital, maxillary, 
orbital and zygomaxillary breadths: JUB, FMB, EKB, MAB, OBB, ZM B; basion-nasion 
length: BNL) we observe small and gracile (left) to large and robust (right) morphologies, 
while vertically along the axis of PC 2 (positively correlated with parietal and frontal 
breadth: XPB, XFB; negatively correlated with position of minimum frontal breadth: 
WFBZFS, palatal length: PALL, dental size: BLM 1 and BLM2, and robusticity: ST, IN, 
TR, RO, SK) the broad crania (upper half) are separated from the relatively narrow skulls 
with large teeth (lower half). Three interesting points are disclosed by the geographical 
distribution seen in Figure 11. First, the group of crania in the upper left quadrant that may 
be described as broad, short and gracile, is mainly composed of Asians (Southeast and East 
Asians 52-8%), Europeans (16 7%), and to a lesser degree, Africans (13-9%). Second, the 
three skulls furthest to the left along PCI are an Andaman Islander, a bushman and an 
individual from Sarawak. They show a very extreme morphology in relation to all other 
crania, and through their diverse origins, show that extreme gracility in terms of both size 
and shape, has been achieved in several modern human groups. And third, it is possible to 
identify the broad-robust (right upper quadrant) from the narrow-robust crania (right lower 
quadrant). The combination of large size, breadth and robusticity may be found in 
individuals from all populations except the Australians, although the majority is formed by 
the Afalou/Taforalt (227%), many of the Fueguian/ Patagonian (2(M%), and many African 
(13-6%), European (13-6%) and East Asian skulls (13-6%). A few skulls from diverse 
regional origins (two Southeast Asians, two African, one Fueguian, one Polynesian and the 
Zkn UC 103 fossil), are long and robust but comparatively narrow, although the differences 
with the group of long, broad and robust crania are not very pronounced. H owever, all the 
Australian crania are considerably narrower than the rest of skulls analysed, forming a clear 
cluster at the bottom of Figure 11. Australian crania depart from the main pattern in that 
pronounced robusticity is observed not only in very large skulls. The group of crania 
identified as G in Figure 7 is not particularly large, and yet presents high levels of 
robusticity. Their distinctive feature is the small width of the vault. 

In summary, the results of the previous analyses show that: 

(1) There is a close relationship between overall cranial robusticity (including facial and 
occipital traits) and cranial size whereby variation in the metrical dimensions has an effect 
on the expression of robusticity. 

(2) There is an even closer relationship between facial robusticity and craniodental 
dimensions. 

(3) T he main metrical factor underlying these relationships is size, in particular cranial length, 
upper facial breadth, facial- basal length, palatal size and dental dimensions. 

(4) A secondary metrical factor is related to cranial shape, contrasting broad gracile skulls with 
narrow robust ones. 

(5) T he integration of the two components of variation, size and shape, disclose 12 anatomical 
combinations with differing levels of robusticity. 
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(6) There is a clear geographical pattern to craniofacial robusticity among recent humans, 
whereby the robust skulls from Australia, Tierra del Fuego, Patagonia, Afalou and Taforalt 
are contrasted with the gracile character of most recent populations. 

(7) M ost robust individuals present absolute large size, in both length and breadth of the vault, 
except for Australian crania that are robust and narrow. 



The results obtained in this study indicate that spatial relationships of craniofacial form are the 
most important correlates to robusticity. The amount of variance in robusticity (taken as a 
single complex) explained by metrical dimensions varied from 25-40%. Two anatomical 
relationships of robusticity are identified, one in terms of cranial size and another one in terms 
of cranial shape. The relationship with the group of variables that represent size is the stronger 
of the two. The anatomical parameters of the complex representing cranial size remained 
stable in all the analyses performed, and may be characterized as comprising cranial length, 
upper facial breadth, facial-basal length, palatodental dimensionsand basal breadth. The level 
of cranial and facial robusticity (expressed in terms of the size of the supraorbital ridges, 
zygomatic trigone, zygomaxillary tuberosity, occipital torus, depth of the infra-glabellar notch 
and sagittal keeling) increases as the variables within the size complex increase in value. T his 
relationship identifies a group of skulls that are either very robust or very gracile (upper and 
lower quadrants of Figures 3 and 5), the robust crania being on average 8-12% larger than the 
gracile ones. T hese two groups of skulls, very robust and very gracile, do not differ in relation 
to mean cranial and frontal breadth, although all gracile crania are relatively broad (XPB: 
CV =44), whereas robust skulls may be broad or narrow (X PB: CV =7-3). The relationship of 
cranial robusticity with a shape complex, identified through the PCA, contrasts gracile crania 
with broad frontals and parietals with robust narrow skulls that also have very large teeth, 
palates and a relative backward position of minimum frontal breadth. This relationship is 
independent of size as described above. 

The combination of the size and shape factors that constrain the expression of cranial 
robusticity define a number of anatomical groups of skulls that follow specific parameters in 
terms of cranial breadth, size and robusticity. Two different anatomical circumstances in 
which cranial architecture needs superstructures can be recognized (a) in very long skulls with 
broad upper faces and large palates (independent of the breadth of the vault), and (b) in very 
narrow skulls with large teeth (independent of cranial length and upper facial breadth). The 
group that presents the greatest robusticity is one that combines these two conditions, i.e. very 
long and narrow skulls, with broad upper faces and large teeth. It is possible to classify crania 
into robust and gracile groups, which may be divided into six categories: 



Discussion 



Functionality 



Gracile crania 
very small 

small and very narrow 
small and narrow 
small and ± broad 
small and broad 
very broad 



(a) very large 

very large and very broad 
large and broad 
large and ± narrow 
very large and very narrow 



(b) very narrow 



Robust crania 
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Studies on the relationship of each feature of robusticity and craniodental dimensions 
through discriminant analyses found that the expression of each trait may be predicted with 
ca. 100% accuracy (Lahr, 1992, 1996). The difference with the results obtained here, which 
reveal a relationship common to all cranial superstructures, suggest that the functional 
size/ shape constraints are strong, but vary in their determination for each trait. Although this 
study did not explore the biomechanical relations between metrical dimensions and the 
superstructures, it shows through the redundancy values that it is the dimensions of a skull that 
affect the expression of robusticity, while the latter has little effect on the metrical parameters. 
T he close relationship between palatodental size and facial robusticity is clearly seen in the plot 
of variables along the first and second principal components of variation (Figure 6). This 
association gives support to the functional interpretation of the facial superstructures in terms 
of resistance to masticatory stress and the potential effects of localized stress on the 
development of a particular superstructure (E ndo, 1970; R ussell, 1983, 1985; D ernes & C reel, 
1988; Hilloowala & Trent, 1988a,b). Russell (1985) argued that the moment arm of the 
downward pulling force of the temporalis is represented by the distance from mid-face to the 
temporal fossae. Russell (1985) used minimum frontal breadth to represent this dimension. It 
is possible that bifrontal breadth, biorbital breadth and bijugal breadth represent equivalent 
dimensions, by which the breadth of the upper face relates to the necessary strength to produce 
effective muscle forces, and affects robusticity as part of the size complex identified in this 
study. 

T herefore, the results from the present study could be interpreted as suggesting that each of 
the features of robusticity that form the cranial superstructures has a specific relationship with 
the muscular and masticatory forces acting on the skull, but that they also form a single 
functional complex, with common underlying anatomical constraints that explain about 40% 
of the variation observed. 

Regionality 

Among modern populations, there is a marked regional pattern of expression of cranial 
robusticity and its anatomical correlates— cranial length and breadth, upper facial breadth and 
palatodental size. Although Figures 3 and 5 show that recent modern variation in robusticity 
forms a gradient from extremely gracile to extremely robust, the distribution of populations 
along this gradient (Figure 9) is not gradual. T his becomes clear when the percentage of crania 
from each region in each of the four quadrants of Figures 3 and 5 is calculated (Table 4). The 
results confirm visual observations of robusticity, and separate recent populations on the basis 
of integrated metrical and non-metrical parameters: the large and robust groups— A ustralo- 
Melanesians and Fueguians-Patagonians— from the small and gracile ones— sub-Saharan 
Africans, Southeast Asians and EastAsians. The Europeans fall in between the extremes, with 
a relatively large number of cases of robust facial features, although not presenting any case of 
pronounced robusticity (i.e. circled in the upper right quadrant of Figures 3 and 5). The 
sample sizes of Polynesians, Inuits and Ainu are too small to allow characterizations of these 
groups. H owever, the findings again confirm visual observations, with the New Zealand and 
Ainu crania falling within the parameters of the robust populations (upper right quadrant), and 
Inuit within the gracile ones (lower left quadrant). Of the Pleistocene and early Holocene 
modern fossils included, the majority of the Afalou/ T aforalt are clearly robust, as are the 
Zhoukoudian UC 101 and 103 remains, whereas the M iddle Eastern late Upper Palaeolithic 
(O halo, Nahal Ein Gev), and especially the few complete N atufians, tend to be more gracile. 
A further distinction can be made. As mentioned above, cranial robusticity is expressed in very 
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large skulls (a) and in very narrow ones (b). In terms of regionality, the striking feature is 
that the separation of modern crania along the large-robust/ small-gracile dimensions is 
applicable to all populations but the Australo-M elanesians. The latter population follows 
certain parameters of its own, by which robusticity becomes relatively decoupled from cranial 
size to follow a closer association with cranial breadth and dental size. In this respect, 
Australo-M elanesian crania are different from most modern human skulls. 

Although sexual dimorphism in most modern groups is expressed as different levels of 
robusticity in males and females, the different ranges of variation in robusticity observed in two 
series of sexed skulls shows that within-population variation due to sexual dimorphism is less 
important than the general metrical and non-metrical parameters of the population as a 
whole. 

P hylogenetic interpretations 

Although ontogenetic bone remodelling due to localized biomechanical stress probably affects 
the full development of cranial superstructures, the distinct regional patterns in robusticity 
cannot be explained solely by differential dental loading or para- masticatory use of the 
dentition in currently robust populations. I n his comment on R ussell's work, T urner 1 1 (1985) 
argued this point clearly by stating that the populations that load their teeth the most— the 
Eskimos— are not amongst the more robust. Furthermore, recent studies on in vivo bone strain 
analysis have shown that masticatory forces do not create sufficient strains on the supraorbital 
area of primates to account for supraorbital torus development (Hylander et al., 1991). 
M univariate analyses of modern human variation have shown that main cranial dimensions 
reflect the genetic differentiation of modern groups throughout the Upper Pleistocene. 
Therefore, it may be assumed that variation in these dimensions, which would affect the 
expression of cranial robusticity, are largely of a genetic basis. It is suggested that cranial size 
and robusticity, and to a lesser degree cranial shape and robusticity, form functional complexes 
that co-varied along the evolution of present cranial diversity through varying dimensional 
changes in the different populations. 

Accepting that size/ robusticity form an evolutionary unit of covariance, the present 
regional pattern of robusticity could still be the result of either phylogenetic or purely 
functional processes. Functionally, it is possible that as dental size and the size of the 
muscles of mastication (especially the size of the temporalis which is closely associated to 
cranial length) reduced during the late Pleistocene, an entire functional complex was 
released, and robusticity was lost in all populations that underwent this process. If this was 
the case, then the distribution of robusticity and its anatomical correlate of size does not 
carry any phylogenetic information, because it could have occurred in parallel in most 
modern populations. However, the association of loss of robusticity with technological 
developments is not straightforward. Among hunter- gatherer groups like IKung bushmen, 
Inuit, Australians, Onas or Tehuelches, the whole range of levels of robusticity are 
observed. Brace (1979) suggested that cooking was the main technological development that 
released the pressure on the masticatory apparatus, but the use of fire and some 
pre-masticatory preparation of food apply to all recent populations, gracile or robust. 
Furthermore, dental reduction in southeastern Australia during the early Holocene occurs 
without any indication of technological developments (Brown, 1992). Therefore, the 
incongruence between the mechanical sources of stress in a population and its expression of 
robusticity would suggest that mechanisms other than the release of a selective pressure 
acted upon the gracilization of the modern skull. 
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If craniodental dimensions are in fact more genetically constrained, its reduction in later 
Pleistocene populations would have resulted from specific adaptations, which although 
homologous in most regions because of similar physical constraints of the skull and 
technological developments, would be indicative of phylogenetic relationships among modern 
groups through their spatial and temporal pattern. This is consistent with the results of 
multivariate analyses, which show that although most recent populations have undergone a 
size reduction in relation to earlier modern fossils (Kidder et al., 1992; Lahr, 1996), a unique 
combination of cranial dimensions allows the identification of most human groups (H owells, 
1973, 1989; Brace et al., 1990; Pietrusewsky, 1990, 1992; Wright, 1992). 

If an epigenetic basis for cranial superstructures is accepted, cranial gracilization in most 
recent populations and on the other hand, the retention of robusticity in a few groups, could 
be examined through the fossil record of modern humans. 

Although the definition of modern cranial form is a controversial subject (Day & Stringer, 
1982; Wolpoff, 1986; Brown, 1990; Lahr, 1996), it is generally agreed that a modern skull 
should present a relatively small face tucked under the vault that is relatively short and high, 
a relatively vertical forehead, parietal enlargement, a relatively rounded occiput, a flexed 
cranial base, a canine fossa, an occipital protuberance in the occipital bone and mental 
eminence or chin. M ost descriptions of modern H . sapiens would also include skeletal gracility 
as characterizing the group. The characterizations described above relate mainly to cranial 
shape rather than size, and are observed in a group of skulls from the early U pper Pleistocene 
of Africa and Israel at the exclusion of other hominid populations of possible similar age, like 
the Neanderthals or the Ngandong sample (Stringer, 1993). However, these early modern 
fossils are considerably larger and more robust than most recent humans. K idder et al. (1992) 
estimate that the early modern sample of Skhul and Q afzeh and the early U pper Palaeolithic 
hominids of Europe are from 10-30% larger in main cranial dimensions than recent 
Europeans. The larger size and more pronounced robusticity in the early modern sample are 
consistent with the results obtained in this study, by which the metrical dimensions of all early 
modern fossils imply higher levels of cranial robusticity than most recent populations. The 
ancestral condition for early modern humans is one of large size (in both breadth and length) 
and robusticity, and the later differentiation of populations must have occurred partly through 
reduction along these parameters. It follows that a modern cranial "shape" appeared before 
cranial size reduction and gracilization occurred, and a definition of a modern human skull 
cannot include cranial gracility among its parameters. 

There is still the case of the distinctiveness of the Australian sample. Table 5 presents the 
value of the metrical variables with moderate to high loadings on the first and second principal 
components in a group of late M iddle Pleistocene and early modern fossils. When comparing 
the values in Table 5 with those of Table 3 (the dimensions of the 12 anatomical groupings 
defined by the plot of the first and second PC s), it is evident that early modern fossils are large 
and robust, but also broad. Although they are not complete enough to be included in the 
multivariate analyses, they present values which are most consistent with the groups in the 
upper right quadrant of Figures 7, 8 and 11 (groups L or K). Therefore, the distinctive 
anatomical combination present in Australian crania of a very narrow vault and pronounced 
robusticity, does not represent a plesiomorphic state. They do not share these anatomical 
parameters with either the early modern sample from Africa/ Israel nor with the Ngandong 
sample from Java. C urrent dates for the occupation of Australia suggest a very early entry into 
the continent (Roberts et al., 1990), and both archaeological and morphological data are 
compatible with a subsequent period of relative isolation (K irk, 1981). A possible explanation 
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for the Australian morphology is that during a comparatively long period of isolation and 
varying demographic parametres (Smith & Sharp, 1993), this population underwent significant 
morphological differentiation as it adapted to late Pleistocene Australia. In terms of 
evolutionary pathways, diversification subsequent to the original adaptations would be 
responsible for the complex mosaic pattern of constancy in some features and pronounced 
heterogeneity in others observed among more recent aboriginal groups. 

Implications for modern human origins 

Because several of the features of robusticity examined in this study have been advocated as 
evidence for regional continuity in southeast Asia, the results obtained have important 
implications for the on-going debate about the origins of modern humans. 

First, it was shown that the cranial superstructures respond to a functional complex, i.e. 
although they may vary in their total independent expression, their incidence is correlated and 
they respond to common underlying factors. This means that, on a worldwide basis, a skull 
with large supraorbital ridges is likely to have large zygomaxillary tuberosities, large zygomatic 
trigones, deep infraglabellar notches, and even occipital tori. These results strongly reject 
the use of each of these features independently as phylogenetic markers (Wolpoff et al., 
1984; Frayer et al., 1993). Furthermore, these results also invalidate suggestions that the 
combined occurrence of all these markers in Australian skulls carries phylogenetic weight 
(H abgood,1989), for it is now shown that the expression of these features is correlated. 

Second, the geographical pattern of cranial robusticity among recent humans shows that 
although most populations have undergone a process of reduction in cranial size accompanied 
by loss of cranial superstructures, this process has not been universal. Furthermore, Pleistocene 
modern fossils are comparatively larger and more robust than most recent humans, a pattern 
that is consistent with the relationship between cranial size and robusticity identified in this 
study. Therefore, both the spatial and temporal patterns of size and robusticity highlight the 
need to disassociate the definition of a modern cranial form from levels of robusticity. 
Arguments that early U pper Pleistocene African fossils are not "modern" based on high levels 
of robusticity do not apply. 

And last, these analyses show that Australian cranial morphology, at the centre of the debate 
about modern human origins, is indeed an outlier among modern anatomical patterns. 
H owever, the results also suggest that the Australian pattern does not result from the retention 
of a plesiomorphic combination of cranial dimensions and robusticity that would associate 
them to the Ngandong hominids. Although the cranial superstructures are in themselves 
plesiomorphous for all modern humans, the anatomical parameters that constrain their 
expression in Australian crania are a characteristic of this group. The most likely explanation 
for their distinctiveness probably lies in the early occupation of Australia and subsequent 
differentiation within the continent. Wright (1976) suggested a similar evolutionary mechanism 
of local response to selection once in Australia to explain large dental size in the late 
Pleistocene Aboriginal population of this continent. 

Conclusions 

The present study was carried out to examine the existence of a relationship between cranial 
metrical dimensions and levels of robusticity. The results obtained show that many of the 
superstructures of the skull have a correlated expression, so that cranial robusticity may be 
treated as a morphological complex by which the probability of a skull presenting not one but 
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several superstructures is high. The canonical correlation analyses showed that there is a 
relationship between this morphological complex "robusti city" and cranial metrical dimensions, 
which remained constant as the major component in all analyses. Furthermore, the analyses 
show that although the dimensions of the skull affect the development of cranial superstruc- 
tures, the latter have a much less important effect on the metrical parameters. The main 
correlate to robusticity is cranial size, in particular cranial length, upper facial breadth, 
facial-basal length and palatodental dimensions. The results from the canonical correlations 
allow the identification of robust and gracile groups of crania, the former between 8-12% 
larger than the latter. PCA revealed that there is also a less important relationship between 
cranial robusticity and shape of the skull, contrasting broad gracile skulls with narrow robust 
ones with large teeth. An orthogonal plot of the cranial sample along the first two principal 
components shows that the combination of the size and shape parameters separates the broad 
from the narrow robust skulls. 

A canonical correlation analysis using only the features of facial robusticity obtained a 
stronger relationship with size than when including all the superstructures studied. These 
results, together with previous findings of stronger relationships between each cranial 
superstructure and metrical parameters, suggest that although there is a "robusticity" complex 
with common underlying parameters, each of the features within this complex responds to 
specific cranial dimensions and mechanical forces for their full expression. Although both 
phylogenetic and ontogenetic processes act through developmental mechanisms based on 
functional relationships, the correlation of robusticity with dimensions that characterize 
different human populations suggests that there are strong phylogenetic processes involved in 
determining the expression of cranial superstructures. 

Consistent with the results obtained from recent modern humans, early modern fossils show 
a combination of great robusticity and large size along the parameters identified in these 
analyses, both in terms of cranial length/ upper facial breadth and cranial breadth. This 
combination is still present in relatively late Pleistocene fossils, like the Zhoukoudian Upper 
Cave and the Afalou/ T aforalt samples. Among recent modern populations, most groups have 
undergone a process of gracilization in size and superstructures, with a few exceptions like 
Fueguian, Patagonian, some Ainu and to a lesser degree some European crania, which still 
show great size and relative robusticity. Australian crania are also robust, but differently from 
the above mentioned groups, robusticity in this group is associated with very narrow skulls, 
some of which are not particularly large. This distinctive pattern in Australia is not related to 
the retention of an ancestral condition because it is not found in conjectural ancestral forms, 
whether the early modern African/ M iddle Eastern samples or the Ngandong remains. 

T he correlation between robusticity and cranial dimensions highlights the need to study the 
cranial superstructures within the context of metrical variation, and not individually as 
evidence of phylogenetic relationships. 
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